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I. 


INTRODUCTION 


Electromagnetic  scattering  from  a moving  body  generally  is  a very 
difficult  problem  to  solve.  Due  to  the  motion  of  the  scatterer,  the  con- 
stitutive parameters,  y,  e and  a are  not  only  functions  of  position  but 
also  functions  of  time,  thus  introducing  modulation  in  the  scattered 
fields.  To  make  things  worse,  the  boundary  conditions  for  a moving  body 
are  much  more  complicated  than  those  encountered  in  problems  involving 
stationary  bodies  only  [1].  These  are  probably  the  main  roadblocks  to 
obtaining  the  solutions  of  the  scattering  problems  involving  moving  bodies. 

However,  in  some  cases,  approximate  solutions  of  the  scattering 
problems  involving  moving  bodies  can  be  obtained.  In  this  report  we 
shall  discuss  the  backscatter  of  a large  rotating  conducting  cylinder. 
Approximate  solutions  of  the  backscatter  will  be  obtained  using  the  quasi- 
stationary method.  The  investigation  of  this  problem  stems  from  the  cur- 
rent interest  in  the  DADAM  phenomena  which  involve  the  identification  of 
moving  vehicles  based  on  the  characteristic  modulations  that  they  impose 
on  their  radar  returns. 


II.  FORMULATION  AND  NUMERICAL  SOLUTIONS 

An  infinitely  long  conducting  cylinder  is  rotating  at  an  angular 
velocity  a about  its  axis.  A plane  wave  is  incident  normally  upon  the 
cylinder.  Due  to  the  rotation  of  the  cylinder,  the  constitutive  para- 
meters are  then  periodic  functions  of  time  with  period  T.  For  a cylinder 
with  arbitrary  cross-section,  T=2Ti/a.  However,  for  an  orthogonal  polygon- 
cylinder,  T=27r/Na,  where  N is  the  number  of  sides  of  the  polygon.  Assume 
that  1/T  is  much  smaller  than  the  frequency  of  the  incident  wave.  Then 
the  backscatter  is  both  amplitude  and  phase  modulated  which  can  be  written 
as 


E 


s 


L{t) 


e 


(1) 


where  wq  is  the  incident  angular  frequency  and  5:(t)  is  a complex  function 
of  t with  period  T.  Because  of  the  assumption,  2t^/T<<wq,  the  negative 
frequency  components  of  the  Fourier  transform  of  E^  will  be  negligible. 
Therefore  E^  can  be  considered  as  an  analytic  signal. 

An  exact  solution  of  z(t)  can  be  obtained  in  two  ways.  The  first 
is  to  solve  the  scattering  problem  in  the  time  domain  using  boundary 
conditions  for  moving  bodies.  The  second  involves  transformations  from 
the  laboratory  frame  to  a non-inertial  frame  which  is  corotating  with 
the  cylinder.  Both  methods  are  very  tedious  and  solutions  difficult  to 
obtain.  However,  under  the  assumptions  that  WQ>>2Tr/T  and  that  the  linear 
velocity  of  the  cylinder  (v  = aa,  a = the  largest  radius)  is  much  smaller 
than  the  velocity  of  light  in  vacuum,  an  approximate  solution  of  r(t) 
can  be  much  easily  obtained  using  the  quasi-stationary  method.  Using  this 


1 


method,  the  backscatter  i;(t)  at  each  time  t is  obtained  as  if  the  cylin- 
der were  stationary.  The  ordinary  frequency  domain  techniques  and  the 
boundary  conditions  for  stationary  bodies  can  then  be  applied  to  solve 
the  problem.  The  inequalities  uq»2it/T  and  c>>aa  are  easily  satisfied 
with  high  incident  frequency  and  low  rotating  velocity.  Thus  the  quasi- 
stationary method  can  be  employed  in  many  practical  situations. 

In  applying  the  quasi-stationary  method,  GTD  (geometrical  theory 
of  diffraction)  techniques  can  be  employed  for  solving  i(t)  from  a large 
conducting  cylinder.  In  the  following,  we  shall  first  briefly  review 
the  GTD  techniques  pertinent  to  this  problem.  For  a more  detailed  form- 
ulation readers  are  referred  to  a paper  by  Wang*  [2]. 

As  shown  in  Figure  1,  a conducting  convex  cylinder  is  illuminated 
by  a line  source.  The  cylinder  is  modeled  as  a polygonal  cylinder  with 
N flat  faces.  The  line  source  illuminates  the  polygonal  cylinder  giving 

LINE  SOURCE 


Figure  1.  The  self-consistent  representation  of  the 
surface  waves  traveling  on  the  triangular 
cyl inder. 


*The  authors  wish  to  thank  Dr.  N.  Wang  for 
program  and  the  reproduction  of  Figures  1 


the  use  of  his  computer 
and  2 in  this  report. 
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rise  to  reflected  fields  and  edge-diffracted  fields,  which  consist  of 
singly-  and  multiply-diffracted  fields.  The  multiply-diffracted  rays 
are  due  to  those  surface  waves  traveling,  at  least  once,  along  the  sur- 
face of  the  polygonal  cylinder  and  then  diffracting  from  the  edges  of 
the  polygon.  The  surface  waves  on  each  side  of  the  polygon  are  conveni- 
ently combined  into  two  such  waves  traveling  in  opposite  directions. 

Hence  there  are  totally  2N  surface  waves  on  a polygonal  cylinder  with  N 
sides.  The  amplitudes  of  these  surface  waves  need  be  solved  first.  We 
shall  use  the  triangular  cylinder  of  Figure  1 as  an  example  to  set  up 
equations  for  solving  these  amplitudes. 

In  Figure  1 the  C-j  are  the  unknown  amplitudes  of  the  surface  waves 
on  the  polygon  faces.  The  surface  wave  field  is  defined  as 

p-jkp 

U{p)  = Ci  ^7^-  , (2) 

/p 

where  k = ^ is  the  distance  from  the  caustic  of  the  surface 

wave  to  the  field  point.  This  field  equation  is  not  valid  at  the  caustic. 
The  boundary  condition  at  the  corners  may  be  expressed  as  follows: 


at  corner  A, 

^1  = C3TCA+C4RBA+V1  . (3) 

^6  = CAi'BA'^CgRcA+Ve  . (4) 

at  corner  B, 

^2  = CiTaB'^C5'^CB'^'^2  ' 

^4  = C5'’’CB'*’^1'^AB‘^V4  , (6) 

at  corner  C, 

C3  = C2TBC+C6RAC+V3  . (7) 

C5  = C6TAC+C2RBC+V5  • 


Equations  (3)  - (8)  can  be  written  in  matrix  form, 

[Z](C)  = (V)  , (9) 

where  [Z]  is  the  coupling  matrix  and  (V)  is  the  excitation  column. 

The  coupling  matrix  [Z]  is  unsymmetric.  Its  elements  specify  the 
interactions  among  the  surface  waves.  The  diagonal  elements  are  Z-ji  = l. 
The  other  elements  are  0,  T^y  or  Rxy.  The  transmission  coefficient  T^y 
and  the  reflection  coefficient  R^y  af'e  given  by 
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(10) 


-jkp, 


xy 


J Y \J  1 

T„..  = ® ^ 1 D3(360°-ay,0°.90°) 


xy 


-jkp 


\y  = i D^(0°.0°,90°) 


xy 


2 S' 

h 


(11) 


where  p^^  is  the  distance  between  the  corners  X and  Y,  a„  is  the  wedge 
angle  for  corner  Y and  D^(<t>.4'' .6q)  is  the  scalar  wedge  diffraction  coef- 
ficient given  by  Kouyoumjian  and  Pathak  [3]. 


The  elements  of  the  excitation  column  (V)  relate  the  source  excita- 
tions to  the  surface  waves  on  the  polygon  surfaces.  The  typical  element 
of  (V)  is  given  by 


Ds(0°.<^’,90°) 

h 


(12) 


where  p^^^  is  the  distance  between  the  source  and  corner  X.  However, 
when  a corner  is  not  directly  illuminated  by  the  source,  the  excitation 
element  associated  with  that  corner  is  0.  Furthermore,  for  the  TM  case, 
the  diffraction  coefficient  Ds  in  Equation  (12)  is  identically  zero  and 
should  be  replaced  by  its  normal  derivative  [2]. 


From  the  matrix  equation.  Equation  (9),  we  obtain  the  amplitudes 
of  the  surface  waves  on  the  polygon  surfaces.  The  scattered  field  at  an 
observation  point  P is  equal  to  the  sum  of  the  fields  associated  with 
all  the  rays  emanating  from  the  cylinder  and  passing  through  P: 


U(P)  = I Un(P) 
rays 


(13) 


which  includes  the  reflected  field  and  the  diffracted  field.  Referring 
to  Figure  2,  the  reflected  field  associated  with  the  reflected  ray  S'P 
emanating  from  the  image  S'  is  given  as 


£jkp^ 

*^^5'  P 


(14) 


where  the  positive  sign  is  for  the  Neumann  (hard)  boundary  condition 
associated  with  the  magnetic  line  source,  and  the  negative  sign  is  for 
the  Dirichlet  (soft)  boundary  condition  associated  with  the  electric 
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c 

IT 


Figure  2.  Direct  ray  5P,  ref.iected  ray  S'P,  singly- 
diffracted  ray  SAP,  and  surface-wave 
diffracted  ray  ABP. 


line  source.  The  singly-diffracted  field  associated  with  the  ray  SAP 
is  given  by 

. -jkPSA 

• (15) 

•^PSA  h ^AP 

and  the  surface-wave-diffracted  field  associated  the  surface  wave  C]  is 
given  by 

. ■j'^PAB  1 p’J^PbP 

•^PAB  h '^Pgp 

These  fields.  Equations  (14),  (15)  and  (16),  exist  only  in  their  illumi- 
nated regions.  In  the  shadowed  regions  the  corresponding  fields  vanish. 
Summing  the  reflected  fields,  the  singly-diffracted  fields  and  the  sur- 
face-wave-diffracted fields,  we  then  obtain  the  scattered  field  at  P. 
This  field  depends  on  the  orientation  of  the  polygonal  cylinder.  As  the 
cylinder  rotates,  the  scattered  field  at  P,  i.e.,  U(P),  then  varies  with 
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time.  Therefore  one  may  write  U(P)  as  U(P,t)  which  is  essentially  equi- 
valent to  i;(t)  of  Eauation  (1).  A computer  program  is  written  for  calcu- 
lating U(P,t)  (or  i:(t)). 

In  this  investigation  we  are  interested  in  studying  the  power  den- 
sity spectrum  of  the  backscattered  field.  The  power  den^sity  spectrum  of 
the  backscattered  field  can  be  easily  calculated  after  z(t)  of  Equation 
(1)  as  a function  of  t is  obtained. 

Because  z(t)  is  a periodic  function  of  t with  period  T,  it  can  be 
represented  as  a summation  of  its  Fourier  series, 


2(t)  = I F(m)  . (17) 


If  |F(m)|  % 0,  for  |m|^NQ/2,  then  Equation  (17)  can  be  adequately  approx- 
imated by 

Nn/2-1 

r(t)  = ° I F(m) 

m=-No/2 

= F(m-No/2)  . (18) 

m=0 

Let 

jT^Not/T  _ 

lit)  e =Ei(t)  , (19) 


and 


F(m-No/2)  = Fi(m) 


(20) 


Then 


5;^(t) 


No-1 

m=0 


F^  (m) 


j2Trmt/T 

e 


The  inverse  transform  is 


Ft  (m) 


Nq-I 

- 

0 tNg/T^O 


L^(t) 


(21) 


(22) 
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Equation  (22)  can  be  solved  using  FFT  programs. 
From  Equations  (1)  and  (17),  we  have 


Es  = I F(m)  e 


j(wQ+2  m/T)t 


m=-oo 


The  power  density  spectrum  of  Equation  (23)  is  defined  by  [4] 

rT„ 


S(w)  = 

^'o 


-To  m= 


^ _ j[(w-u)„)-2TTm/T]t  1 2 

\ F(m)  e dt 


(23) 


P ” sin[(u-u  )-2TTm/T]T 

= f-  I F(m)  - ° ° 

0 m=-«> 


(u)-Wo)-27Tm/T 


_ sin[(u)-a)o)-2iTm7T]To 


As  To- 

S(u)  = 2-n  |F(m)|  6[(w-w^)-2-irm/T] 

m=-oo 


(24) 


(25) 


Thus  the  power  density  spectrum  of  the  scattered  field  from  a rotating 
cylinder  consists  of  spectral  lines  located  at 


u)  = (i)o'*'2'iTm/T  = u)o"HnNa  , m-0,  + 1 ,±2, • • • • . (26) 

Figures  3 to  10  show  some  power  density  spectra  of  the  backscattered 
fields  from  a large  rotating  conducting  cylinder.  The  cross-section  of 
the  cylinder  is  an  N-sided  orthogonal  polygon.  The  spectra  in  Figures  3 
and  4 are  set  to  0 dB  at  m=0  (u=uo).  The  dB  values  of  the  spectra  in 
other  figures  are  then  adjusted  relative  to  these  two  spectral  lines. 
Examining  these  figures,  we  note  several  interesting  features; 

(1)  The  spectrum  is  symmetric  about  m=0.  This  property  is  due 
to  the  symmetry  possessed  by  the  scattering  orthogonal 
polygon-cylinder.  In  the  next  section,  we  shall  show  that 
the  power  density  spectrum  of  the  backscattered  field  from 
a rotating  cylinder  with  an  arbitrary  cross-section  in 
general  does  not  have  this  symmetry  property. 
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Figure  4.  Power  density  spectrum. 
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Figure  6.  Power  density  spectrum. 
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Figure  8.  Power  density  spectrum. 
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(2)  The  spectrum  at  m=0  is  not  always  a maximum  as  evidenced 

from  Figures  7 and  8.  This  is  probafy  true  when  the  length 

of  the  polygon  side  exceeds  several  times  of  the  incident 

wavelength. 

(3)  As  |m|  increases,  the  spectrum  usually  decreases  monotonical ly 

and  faster  for  the  TM  case.  However,  for  the  TE  case  the 

spectrum  has  a local  maximum  around  |m|  ^ a)r>/Na,  where  wq  is 

the  Doppler  angular  frequency  given  by  ^ 

ojQ  = 2Ti'2aa/XQ  = 2ao(i)Q/c  , (27) 

where  a and  aa  are  tiie  largest  radius  and  the  largest  linear 
velocity  of  the  cylinder,  respectively.  For  example,  in 
Figure  4,  a)Q%20Tra  and  a local  maximum  occurs  at  |m|  = 6 
% oj[)/Na.  The  differences  in  the  shapes  of  the  spectra  in 
the  TM  and  TE  cases  are  probably  due  to  the  fact  that  in  the 
TM  case,  the  induced  current  on  the* surface  of  the  cylinder 
is  perpendicular  to  the  velocity  of  the  cylinder,  while  in 
the  TE  case  the  induced  current  has  a component  parallel  to 
the  velocity. 

(4)  In  all  these  figures,  although  not  shown  explicitly,  the 
spectrum  drops  sharply  when  lm|  > wp/Na.  In  the  next  section, 
we  present  an  analytical  evaluation  of  this  phenomenon. 

Figures  3-10  are  the  backscatter  power  density  spectrum  when  a 
plane  wave  is  incident  upon  a rotating  orthogonal  polygon-cylinder. 

They  can  be  obtained  by  using  a line  current  lO^Xg  away  from  the  center 
of  the  cylinder.  At  such  a large  distance,  the  curvature  of  the  incident 
wavefront  at  the  cylinder  is  very  small  and  the  incident  wave  can  be  con- 
sidered as  a plane  wave.  However,  when  the  line  current  is  only  IOOAq 
or  less  away  from  the  center  of  the  cylinder,  the  line  current  is  in  the 
near-field  zone  of  the  cylinder  and  the  wave  incident  upon  the  cylinder 
is  a cylindrical  wave.  The  monostatic  scattering  power  density  spectra 
when  the  line  current  is  25Xq,  50Xg  and  IOOXq  away  from  the  center  of  the 
cylinder  are  found  to  be  essentially  the  same  as  those  shown  in  Figures  3 
to  10.  Thus  the  power  density  spectra  of  a rotating  cylinder  can  be  ob- 
tained using  near-field  measurements. 


III.  ANALYTICAL  APPROACH 

Numerical  solutions  reveal  some  interesting  properties  of  the  power 
density  spectrum  of  the  backscattered  field  from  a large  rotating  cylin- 
der. In  this  section  we  shall  investigate  this  problem  analytically  to 
further  substantiate  the  results  established  from  numerical  solutions 
presented  in  the  last  section. 

The  vector  potential  due  to  a current  source  J(r,t)  is  given  by 
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(28) 


A(r.t)  = ^ 


5(t'+R/C-t)  dx'dy'dz’df 


where  R = |r-r'l  and  6 is  a Dirac  delta  function.  Taking  the  Fourier 
transform  of  Equation  (28),  ‘ 


A(r.u,)  = ^ 


f[  g-jw(t'+R/c)  dx'dy'dz'dt' 


(29) 


In  the  two  dimensional  problem  considered  here,  J(r',t')  is  independent 
of  z'.  Therefore  performing  the  integration  over  z'  in  Equation  (29), 
we  obtain 


A(p 


1 - — f 

,u))  - 4j  J 


(2) 


J(p',t')  HQ'"'(ioRi/c)e'-^“^'  dx'dy'dt 


(30) 


where  R]  = ip-p'|  and  is  a Hankel  function  of  the  second  kind.  In 
deriving  Equation  (30),  use  has  been  made  of  the  identity 


i; 


-jkR 


4irR 


dz'  = 


(2) 

Ho  (kRi) 

4j 


(31) 


with  R = [R^  + (z-z')2]’/2_  Equation  (30)  becomes 


M e"j’'/A  e~j“p/c  • /^i  ' ~i/-\ 

A(p,u,)  . ° II  J(pM')  e-J"'* 


where  p is  a unit  vector  along  p. 


dx'dy'dt' , 
(32) 


Now  let  J(p,t)  be  the  induced  current  on  the  rotating  cylinder. 
Then  from  the  vector  potential  we  obtain  the  scattered  electric  and  mag- 
netic fields.  In  the  TM  case,  the  scattered  electric  field  has  only  the 
z-component. 

-ia)u„e”3^/A  0-jwp/c  . _ 

" .n„-Ja)(t'-p-p'/c), 


E,(p,a))  = 


J,(p',t')e 


/87ru>p/  C 

And  in  the  TE  case,  the  scattered  magnetic  field  has  only  the  z-component. 


dx'dy'dt' 
(33) 


Ct^OlTlop/C 


[pxJ(p',t')]^ 


dx'dy'dt'  . 

(34) 
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The  factors  outside  the  integrals  in  Equations  (33)  and  (34)  will  be 
designated  as  Tnn  and  T£,  respectively. 

We  shall  evaluate  the  integrations  in  Equations  (33)  and  (34)_by 
transforming  the  coordinate  system  (p,t)  to  the  coordinate  system  (p^.t^) 
which  is  corotating  with  the  cylinder.  The  coordinate  transformation  is 
governed  by  the  following  relations 

X = X]  COSat]  - y-j  sinat]  , 

y = x^  sinat^  + y^  cosat]  , (35) 

t = ti 


These  equations  can  be  derived  by  referring  to  Figure  11.  The  Jacobian 
of  the  transformation  is  unity.  Therefore  Equations  (33)  and  (34)  become 


E2(p,w)  = j J2(p].tj)  e 


-ja)(t-j-p*pj/c) 


dx]dy-jdtj 


_ rr  -jw(tl-'''pl/c) 

H2(p,u))  = T^  [pxJ(pi  .tl)]^  e dx|dy.Jdt^ 


Figure  11.  Transformation  of  coordinates. 
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To  proceed,  we  assume  that  the  incident  plane  wave  has  the  form 
expEjoj^U-x/c)].  which,  if  written  in  the  corotating  frame,  becomes 


exp{ jwQ[t^-(xiCOSat^  - y^sinat^ )/c]} 
= exp{  jwQ[t^-p-|COs(at^+(()l  )/c]} 

.-k 


= e 


I J 

k=-o 


c)  e 


jk(at^+((.^  ) 


(38) 


2 2 1/2  -1 

where  P]  = (x-i  + y-i)  . 4’'|  = tan  (y-i/x-i),  and  J.  is  a Bessel  function. 
In  deriving  Equation  (38),  use  has  been  made  of  tne  relation 

eJxsin^  = ^ j (x)  . (39) 

k=-» 

In  view  of  Equation  (38),  the  induced  current  J(pi,ti)  can  be  written 
as 


- .-k  - . jC(V''“)^i^''^i^ 

d(p]5t'])  " / J ^k  ^1  ^ 

k=-“ 


(40) 


Thus  using  Equations  (36),  (37),  (39)  and  (40),  we  obtain  the  backscat- 
tered  fields. 


E2(p,w)  - T|^ 


“ .-(k+m)  j(k+iTi)i.i 

I I J 3||j^(pi  )Jp^(a)Pl/c)e 

k,m=-” 

-j[u)-uo-(k+m)a]tj 

• e p]dp]d(t)jdtj  , (41) 


H2(p.“)  = 


I I 

k,m=-” 


1_  .-(k+m-1) 


aik(pi)Jm(“Pi/c)e^^''^"'^^i 


f jB^  -j[a)-(OQ-(k+m+l  )a]tj  -jaj  -j[a)-a)Q-(k-Hn-l  )a]tj 


I' 


e e 

p 1 dp  1 d((ii  dt-j 


- e 


(42) 


where  an  and  aj^  are  the  components  of  a parallel^  and  perpendicular  to  the 
z-axis,  respectively.  6|  is  the  angle  between  ai  and  the  positive  Xi- 
axis.  Substitute  m=n-k  into  Equation  (41)  and  m=n-kTl  into  Equation  (42) 
and  carry  out  the  integration  over  t]. 


-:r 
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— / ."fl 

E2(p,oa)  = Zir  Tf^  I I I j 6(u-(0Q-na)  e 
jj 


I 3|i|^(p  ' )Jj^_k(wpi/c)  p-jdp-jd(()-j  , (43) 
k=-“> 


H2(p,u))  = 2ii  J [ - ^ 6(w-a)Q-na)  • J 

n=-"  k=-“ 


j[(n-l  )<}i-j+3(]  j[(n+l  )(j)-i-6i] 

e Jn_l_k(wpi/c)  + e 


p'dp'd*^  . (44) 

The  integrations  in  Equations  (43)  and  (44)  are  over  the  entire 
circumference  of  the  cylinder’s  cross-section.  If  the  cross-section  is 
an  orthogonal  polygon  of  N sides,  then  Equation  (43)  can  be  simplified, 

- , ff  r -n  J"'<’1 

E^(p,u))  = Ztt  T|^  1 1 I j 6(w-u)Q-na)  e e 

n=-“>  i=0 

00 

3„k(p])  Jn_k(“Pi/c)  pjdpjdf]  , 


where  the  integration  is  over  one  side  of  the  polygon  only.  Since 
N-1  jil^ 

I = " - I 


where  6,-  , is  a Kronecker  delta,  therefore  Equation  (45)  reduces  to 

I ) J 

. rr  ? .-mN 


E2(p,(j)  = 2TtN  1^1  I j 6(a)-u)Q-mNa)  e 


I a„k(pi)  Jn,f^_k(<^Pl/c)  pjdpjd*^  . 
k=-<»  (47) 


Similarly,  Equation  (44)  reduces  to 
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oo 


00 


H^Cp.w)  = 2tiN 


. -mN  _ 

I 6{w-w  -mNa)  • I (p') 

(n=-oo  k=-“= 


r j[(niN-l  )(t)'+6^] 

^ j[(mN+l)<t)^-e]]  ] 

+ e '^mV+l -k ^ 1 ^ 7 pjdp-jd((''j  . (48) 


To  solve  Equations  (47)  and  (48),  we  need  to  know  the  quantities 
ail  and  ai.  Unfortunately  these  quantities  cannot  be  easily  found  except 
for  the  case  of  a circular  cylinder.  For  a circular  cylinder  with  radius 
a [5] 


a(ik(Pi ) “ 


-2 


6(p,-a) 


Wl  h(2)(„^p,/c) 


1 


" ^^Pl/c  h(2)’(^^p^/c) 


6(p,-a) 


(49) 

(50) 


which  are  independent  of  (j)] . Therefore  E7(p,w)  and  H2(p,w)  vanish  except 
when  m=0.  This  means  that  the  fields  scattered  from  a rotating  circular 
cylinder  have  the  same  frequency  as  the  incident  wave.  It  is  conjectured 
that  for  an  orthogonal  polygon-cylinder  with  a large  number  of  sides,  the 
a and  ax  will  be  close  to  those  given  in  Equations  (49)  and  (50)  with 
tne  constant  radius  a replaced  by  the  variable  radius  aU-]).  Therefore 
we  shall  substitute  Equations  (49)  and  (50)  into  Equations  (47)  and  (48) 
and  obtain 


^^(p »“) 


= T. 


I N 6(a)-u  -mNa)e  [ —r^ d(t)l  , 

m=-“  k=-»>  H|^  [uoa((t)j  )/c] 

(51) 


H7(p,U))  = Tr 


«>  jmNiJ)-j  ” 1 

^ N 6 (u)-u)Q-mNa)e  ' 


m=-oo 


k=-oc  [wQa(<t.j  )/c] 


2j 


j(6]-<^l) 

-j  (3-j  -<(>■{ ) I 


(52) 


where  we  have  lumped  all  the  factors  nonessential  for  calculating  the 
relative  power  density  spectrum  into  Tm  and  T^.  The  spectrum  consists 
of  spectral  lines.  A single  spectral  line  is  proportional  to  the  square 
of 
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E 


m 


fn/N  jmN(()i 
N e 


“ JmN-kt(‘^o^Na)a(*T)/c] 
L (7\ 

1^=-“  ^[u)Qa(4i^)/c] 


(53) 


or 


H 


m 


,7r/N  jmN((ii 

N e 

-ti/N 


- jsirnt)^ 


J u(2)'r  ^ ! — TTT  ’ 1^°^'*’! 

k— oo  H|^  ^ [a)oa(<))j  )/c]  L 

rvS)a|<^^  )/c  '^mN-kl^(  '^'^1  . 


(54) 

where  we  have  assumed  that  the  range  of  integration  spans  one  side  of 
the  polygon  and  that  3i=ti/2. 


With  low  incident  frequencies,  a)oa(<t>i)/c  <<  1.  The  arguments  of 
the  Bessel  and  Hankel  functions  in  Equations  (53)  and  (54)  are  very  small. 
Thus  the  values  of  E^  and  decrease  very  fast  as  |m[  increases.  This 
means  that  no  significant  spectrum  spreading  occurs  with  low  incident 
frequencies. 

On  the  other  hand,  with  high  incident  frequencies,  ojq  >>  mNa  and 
a(<('i)  > Ap.  The  arguments  of  the  Bessel  and  Hankel  functions  in  Equations 
(53;  and  (54)  are  in  the  order  of  a)[)/2a,  where  wq  is  the  Doppler  angular 
frequency  defined  in  Equation  (27).  If  we  plot  the  Bessel  and  the  inverse 
of  the  Hankel  functions  with  fixed  arguments  against  the  orders,  we  note 
that  if  the  arguments  are  large,  say  > 10,  then  the  values  of  the  Bessel 
and  the  inverse  of  the  Hankel  functions  decrease  very  fast  when  the  orders 
become  greater  than  the  arguments.  This  is  shown  in  Figures  12-15.  Thus 
with  high  incident  frequencies,  we  can  properly  truncate  the  summations 
in  Equations  (53)  and  (54)  at  ±[uQ/a], 


^m 


it/N  jmN(;i'l  JmN-kC“o^('*’l 


Ne 


-it/N 


--[u)Q/a]  hJ^  ^[(*)Qa(<t'1  )/c]  ^ 


(55) 


,ti/N  jmNti 

J _ ...  Ne  I 


-ti/N 


ITT 

k=“[i»)Q/a]  H|^  [i»)Qa(4>1  )/c] 


• |coS(^i  j;N-kC‘-oa('^l)/c]-jsin*i  — a'^'ly/c  '^mN-k^^ V^*l  ‘^^l 
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(56) 


Figure  12.  Plots  of  |Jp(10)|  and  |Jp(50)|. 


It  can  be  easily  shown  from  Equations  (55)  and  (56)  that  in  general 
. This  means  that  the  spectrum  is  generally  not 
A sufficient  condition  for  the  spectrum  to  be  sym- 
metric is  a(((>i  )=a(-<)>i ) for  -u/N  i <)>  i tt/N,  which  can  also  be  derived  from 
Equations  ^5)  and  (56).  For  an  orthogonal  polygon-cylinder,  this  con- 
dition is  certainly  satisfied.  Therefore  the  backscatter  power  density 
spectrum  of  an  orthogonal  polygon-cylinder  is  symmetric  about  m=0. 


^m  ^ ^-ni  ^m  ^ ^-m 
symmetric  about  m=0. 
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Figure  13.  Plots  of  Hn^^(lO)  and  hJ^^(50)  . 

It  is  also  obvious  from  the  properties  of  the  Bessel  and  Hankel 
functions  with  large  arguments  and  large  orders  that  the  values  of 
and  Hp,  are  very  small  for  |m|  > wq/No. 

Table  1 shows  the  comparison  of  the  power  density  spectra  calcu- 
lated using  Equations  (55)  and  (56)  with  those  obtained  in  the  last 
section.  The  agreement  is  very  good.  This  not  only  substantiates  the 
salient  features  of  the  backscatter  power  density  spectrum  of  a large 
rotating  cylinder  established  in  the  last  section,  but  also  indicates 
that  the  approximation  of  using  Equations  (49)  and  (50)  for  polygon- 
cylinders  is  indeed  a good  approximation. 
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Figure  14.  Plots  of  |Jp(10)|  and  |Jp(50)|. 
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TABLE  1 

COMPARISON  OF  THE  POWER  DENSITY  SPECTRA 


* GTD  solutions 
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IV.  CONCLUSIONS 


In  this  report,  two  methods  for  calculating  the  power  density 
spectra  of  the  scattered  fields  from  a large  rotating  conducting  cylin- 
der are  presented.  Results  obtained  using  these  two  methods  agree  very 
well  with  each  other. 

The  backscatter  power  density  spectrum  of  a rotating  cylinder  with 
arbitrary  cross  section  consists  of  discrete  spectral  lines.  For  an 
orthogonal  polygon-cylinder,  the  spectrum  has  the  following  interesting 
features: 

(1)  It  is  symmetric  about  the  incident  angular  frequency. 

(2)  For  the  TM  case,  the  spectrum  decreases  monotonically. 
However,  for  the  TE  case,  it  has  local  maxima  around 
io±uo,  where  wq  is  the  Doppler  angular  frequency. 

(3)  The  spectrum  is  essentially  confined  to  a bandwidth 

wq  -a)[)<u)<a)Q+u[) . 

(4)  The  near-field  spectrum  is  essentially  the  same  as  the 
far-field  spectrum. 
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